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Abstract 

A method for determining the ocean bottom optical albedo, R b , from imwater 
upward and downward irradiance measurements at a shaliow site is presented, tested 
and compared with a more familiar approach that retires additional measurements 
at a nearby deep-water site. Also presented are two new algorithms for the estimatron 

o' ^ fcm ° f the ‘ W “' d "radiance and vertically upward radiance. 

All methods performed well in numencal situations at depths where the influence of 
the bottom on the light field was significant. 
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1. Introduction 


For shallow ocean waters, knowledge of the optical bottom albedo, R b , is required to model 
the underwater 1,2 and above-water 3 light field, to enhance underwater object detection or 
imaging, 8 and to correct for bottom effects in the optical remote sensing of water depth 4,5 
or inherent optical properties (IOP’s). 6,7 Measurements of R b can also help one identify the 
bottom sediment composition, 5 determine the distribution of benthic algal or coral commu- 
nities, 9 and detect objects embedded in the seafloor. Furthermore, values of R b , defined 
as the upward irradiance emerging from the bottom divided by the downward irradiance 
into the bottom, can be used as an integral test for attempted measurements of the bottom 
bidirectional reflectance function. 

Although the value of the irradiance ratio, R(z), equals R b at the bottom, it is not 
possible to measure R(z) right at the bottom, and irradiance measurements just above the 
bottom are difficult to obtain because of instrument self-shadow. An estimate of R b can 
be made by extrapolating to the bottom measurements of R(z) at several depths z near 
the bottom; however, extrapolation is generally unreliable because profiles of R(z) typically 
vary sharply with depth close to the bottom. 10 In vitro R b measurements of small bottom 
samples can be obtained with the method in Ref. 4 or of larger bottom samples with a 
spectral radiometer; however, these are time-consuming processes, the in vitro value of R b 
is not necessarily equal to the in situ value, and it is not clear how representative these 
samples are of larger spatial regions of interest. 

Because in situ estimates of R b from light measurements close to the bottom are subject 
to small-scale horizontal variability of the bottom, it may be preferable to determine R b from 
measurements further from the bottom, thereby obtaining horizontally-averaged values that 
are more appropriate for remote sensing applications and one-dimensional radiative transfer 
modeling. An estimate of R b can be made 9 from a R(z) measurement just below the sea 
surface, R( 0 + ), together with simultaneous measurements in nearby deep water of R(0 + ) 
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and the downward diffuse attenuation coefficient. Similarly, a qualitative algorithm has been 
proposed 5 and tested 11,12 for bottom characterization from remote radiance measurements 
at two wavelength bands simultaneously over both shallow and deep water. Both of these 
methods have the disadvantage, however, that they require that a deep-water site exists 
nearby that has the same water composition, illumination, and surface conditions as the 
shallow- water site. 

A new method of solving the inverse radiative transfer problem for the determination 
of Rb is proposed in Sec. 2. Required are measurements of the upward and downward 
irradiances at one wavelength and at least two mid-water-column depths at only one site. 
Also proposed are two related algorithms for the estimation of Rb from measurements of the 
vertically upward radiance and downward irradiance. These algorithms, like the previously 
developed ones, require knowledge of the measurement distances above the bottom. Results 
of specific numerical tests for all the R b estimation algorithms are presented in Sec. 3, and 
a discussion is given in Sec. 4. 

2. Theory 
A. Preliminaries 


We are interested in the azimuthally averaged radiance L{z y fx) that satisfies, the integro- 
differential radiative transfer equation 

+ cj L(z, n) = b [ ^ (3{z, fi, //) L(z, fi :') d(i' t (1) 

where h and c are the scattering and beam attenuation coefficients, (3 is the azimuthally 
integrated scattering phase function, and /i is the direction cosine with respect to the down- 
ward depth a. All quantities in Eq. (1) implicitly depend on wavelength. The downward 
and upward irradiances are given by 
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( 2 ) 


r\ fO 

E d (z) = 2tc j LiL(z,n)dn and E u (z ) = 2ir J ^ |//| L(z, /j,)dfi. 


The irradiance ratio is 


R(z) = E u {z)/E d {z), 

and R{z b ) = R b for water depth z b . The analogous radiance-irradiance ratio is 


( 3 ) 


R L {z) = irL u (z)/E d (z), 


( 4 ) 


where the vertically upward radiance L u (z ) = L(z, — 1). The factor tt is included in the 
definition of R L (z ) so that R L {z b ) « R b , with R L (z b ) = R b for a Lambertian bottom. 

With increasing depth in optically deep, spatially uniform, source-free waters, R(z), 
R L (z), and the downward diffuse attenuation coefficient, 

1 dE d (z) 


K d {z) = 


( 5 ) 


E d (z) dz 

asymptotically approach values i?oo> Roo, and Kao, respectively. These asymptotic values 
are IOP’s of the water that can be uniquely computed from b, c, and 0. 1 * 


B. Irradiance ratio approach 

A well known model for the irradiance ratio is 4,5,14 

R(z) = Rondiz) + [Rb - R2nd{z )] exp[-2(z b - z)K 2 nd], ( 6 ) 

where the water depth at the shallow site of interest is assumed to be known and the 
subscript 2nd denotes measurements taken at a nearby deep-water site characterized by the 
same illumination, sea surface conditions, and water IOP’s as the site of interest. In the 
derivation of Eq. (6), K is taken to be K(z) = [K d (z) + k { z )}/2 averaged over depth in 
a non-specified manner, where k is the coefficient of attenuation in the upward direction 
of upwelling photons 14,15 (from both in-water scattering and the bottom); a qualitative 
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and quantitative study of k is given in Refs. 4 and 15. In practice, however, the value of 
K 2n d is typically approximated by the vertically averaged K ( [(z) at the deep-water site. 4 A 
rearrangement of Eq. (6) gives the algorithm evaluated by Maritorena et al. 4 for determining 
Rb from irradiance measurements at two sites, 9,16 

Rb = R 2 nd(z) + [R{z) - R 2n d(z )] exp[2(z 6 - z)K 2nd ). (7) 

Equation (7) is often written for 2 = 0 + in hope of applying it to remote sensing applications; 
however, it is valid for any depth, z, and is more accurate at mid-water depths than at the 
surface. A difficulty with implementing Eq. (7) is that a deep-water {2nd) site may not be 
available that matches the water, surface, and illumination conditions of the shallow-water 
site. 

An alternative shallow-water model was previously derived 17 from the radiative transfer 
equation with the eigenfunction expansion method. In this derivation the light field is 
approximated by the sum of two eigenmodes that decrease in magnitude with distance away 
from the surface and bottom, respectively. Expressed in the form of Eq. (7), this model is 

Rb — Roo + [R{z) - Roo] exp [2 ( Zb - z)K 00 ], (8) 

which can be used to estimate Rb from measurements at a single site provided that Zb is 
known and R ^ and can be determined. It is important to note that the asymptotic 
in-water irradiance ratio R^ in Eq. (8) does not represent the same quantity denoted by 
that symbol in Refs. 4 and 9, where R^ is our R 2n d{® + )- Equations (7) and (8) were derived 
with entirely independent approaches; however, their final forms are very similar, and the 
two theoretically converge when R 2n d and K 2n d are measured at large depths in homogenous, 
source-free waters, where R 2nd {z ) = Roo and K 2nd {z ) = A’oo- 

Equation (8) provides a new interpretation of the attenuation coefficient in Eq. (7). The 
derivation of Eq. (7) suggests that K ~ {K d + k)/ 2, whereas the derivation of Eq. (8) 
suggests that K ~ K^. Because n > Kd, 4,15 and therefore K > Kd, there is a question 
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about the appropriateness of taking K to be the vertically-averaged K d (z). Because typically 
k > Aoc > it follows that K = [K d {z) + K.(z)\ 2 nd /2 & K^. Therefore, Eq. (7) may be 
best implemented by taking K equal to the value of K d (z) deep in the euphotic zone, where 
K d (z) « K 0 o, rather than as earlier proposed. 9 ’ 16 This hypothesis is further addressed in 
Sec. 3. 

If Roo and can be determined at the shallow site of interest, then the method of 
Eq. (8) has the distinct advantage over that of Eq. (7) that measurements are required from 
only the one site. One way to determine R ^ and is to calculate them with the procedure 
in Ref. 17 from measurements of b and c obtained, for example, from water samples 18 or with 
a Wetlabs ac-9 instrument 19 . Alternatively, it is possible to estimate R^ and from the 
same irradiance profile measurements (at the shallow-water site) used to form R{z) without 
any direct measurement of the water properties. The value of R^ can be obtained with an 
equation derived 17 from a shallow-water asymptotic approximation of the light field, 

fezM 2 = r6S 

Li + RooJ [E d {z) + E u (z)Y Ijj • W 

To employ Eq. (9), one must subtract and add E u {z) and E d (z ) at two depths, z\ and 
z 2 i square the results, and evaluate the differences between the two depths. Because these 
operations are susceptible to noise, it is important that the irradiance measurements be of 
high quality and their temporal variations be averaged out. The value of in Eq. (8) 
can be estimated as the maximum value attained by K d (z). Because the value of K d {z ) 
is relatively insensitive to Ri,, 1 the value of maxf/CpT)] is typically approximately equal to 
K 0 o. 

C. Radiance-irradiance ratio approach 

If L u (z ) measurements are available rather than E u (z), then /?/, can be estimated from R L (z) 
of Eq. (4) with a new model (derived in Appendix A), analogous to Eq. (8), 


6 


R t = R L + ( R L {z) - R Ll exp(2(z t - z)K, o]. 


( 10 ) 


However, because an equation analogous to Eq. (9) has not been derived for R^, implementa- 
tion of Eq. (10) requires that the value of either be calculated 13 from local measurements 
of the lOP’s or measured in nearby deep water. 

We found with numerical simulations that reasonable estimates of R b can alternatively 
be obtained from R L (z ) with 

R b = R Ld( Z ) + [R L ( Z ) ~ R 2. nd( Z )} ex P[2 (z b ~ z)K 2n d}, (11) 

although we have no analytical justification for Eq. (11) other than its analogy to Eq. (7). 
As with Eq. (7), use of this equation requires measurements at a second (deep-water) site. 

3. Numerical tests 
A. Methods 

Numerical tests were performed to evaluate the accuracies of Eqs. (7)-(ll) for the determi- 
nation of R b . Simulated E u (z), Ed(z), and L u (z) values were generated at 0.25 optical depth 
spacing using the discrete ordinates radiative transfer code DISORT. 20 The surface illumina- 
tion was modeled as a combination of direct collimated sunlight and diffuse skylight. The 
water was defined to have locally homogenous optical properties, a relative index of refrac- 
tion of 1.34 with respect to air, and scattering determined by the Petzold particle scattering 
phase function. 21 Spatially-dependent internal sources, such as from fluorescence, Raman 
scattering, or bioluminescence, were neglected. A Lambertian bottom was assumed, which 
provides a good approximation to the more general, but usually poorly known, bidirectional 
reflectance function. 1,2 Simulations were performed for various values of single scattering 
albedo (u; 0 = b/c), R b , percent direct sunlight, and water optical depth (r b = cz b ). The val- 
ues of R b were taken to be in the range 0 < R b < 0.4, which is consistent with observations 
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in natural waters. 1 

Equation (8) was applied to shallow- water simulations of E u (z) and E d (z) to determine 
bottom albedo estimates, R b (z), as a function of the depth of the corresponding irradiance 
measurements. The values of R <*, and K 0 0 used in Eq. (8) were obtained in two different ways. 
First, /?oo was determined with Eq. (9) and K 0 0 was estimated from either K^z) « Kjz) 
or K 0 o ~ max) [K d (z)\. Second, the IOP’s of the water were assumed to be known from in 
situ measurements, and these were used to calculate R <*, and K^. For comparison, R b (z) 
was also determined from E u (z) and E d (z) with Eq. (7) for combinations of shallow and 
deep-water simulations. Here, K was taken to be, alternatively, K d (z), K d (z) (the vertical 
average of K d (z ) between the surface and the depth of the shallow- water site), and 
[determined from a deep (asymptotic) value of K d (z)\. 

The bottom albedo was also determined from simulations of L u {z ) and E d (z). Equa- 
tion (10) was used to calculate R b (z) from data at only the shallow water site. The values of 
R%o and K 0 0 were calculated from the known water optical properties. In addition, R b (z) was 
determined from L u (z ) and E d (z ) with Eq. (11) for combinations of shallow and deep-water 
sites. The value of K^ n d was determined in the same manner as for the E u —E d approach. 

B. Results 

For all the shallow-water simulations performed, estimates of R b obtained with Eqs. (8) 
and (9) approached the correct value in a nearly linear fashion within the bottom few 
optical depths of the water column. Extrapolation of R b (z ) at two and one optical depths 
above the bottom consistently produced estimates of R b that were accurate to within about 
1%. Given in Table 1 are example R b {z) at three, two, and one optical depths above the 
bottom, as well as the linear extrapolation of the latter two to the bottom. For Table 1, the 
illumination conditions of the simulations were taken to be either overcast (100% diffuse) 
or sunny (75% direct sunlight from a zenith angle of 30°), and in the solutions for R b the 
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values of in Eq. (8) were approximated by ma x[K d {z)\. At a given measurement depth, 
estimates generally improved with increasing value of uq. For example, for R b = 0.2, r b = 5, 
and sunny skies, the error at 3 optical depths above the bottom was 15% for ujq — 0.7 but 
only 3.5% for = 0.9. The value of R b had a small effect on the accuracy of its estimate; 
for large u 0 and small R b , R b (z ) was greater than R b , whereas for small values of u 0 or for a 
combination of large u>o and large R b , R b (z ) was less than R b . Estimates were more accurate, 
but insignificantly so, with overcast conditions. Also, the depth of the water had very little 
effect on the accuracy of R b (z) at a given depth above the bottom; however, in practice, 
instrument noise will be more significant in relatively deep water than in very shallow water. 

In most cases, estimates of R b with Eqs. (8) and (9) were more accurate if in Eq. (8) 
was approximated by ma x[K d {z)\ than if it was replaced by K d {z). This is because large 
values of K 0 0 in Eq. (8) lead to large values of R b , and the estimates of R b were typically less 
than the correct value. For cases where the value of cjq was high and the value of R b small, 
the use of K d (z ) gave slightly, but insignificantly, better results than the use of max[A' ( /(z)]. 

Because Eq. (9) provides only an asymptotic approximation to R^, it was expected 
that estimates of R b would improve if more accurate values of R^, calculated from the 
assumed known water IOP’s, were used in Eq. (8). However, the numerical tests showed the 
reverse to be true; errors in R b introduced by the approximation of Eq. (9) helped counteract 
errors in R b due to the assumptions inherent in Eq. (8). Although Eq, (8) with calculated 
Roc performed similarly in the bottom half of the water column to Eq. (8) with /?<*, from 
Eq. (9), Eq. (8) with calculated /?<*, performed poorly in the top half of the water column 
and even near the bottom slightly underperformed Eqs. (8)-(9). 

For example, shown in Fig. 1 are the estimates of R b as a function of measurement 
optical depth obtained from simulated E u (z) and E d (z) with three different methods: from 
R b (z) = R(z), from Eq. (8) with calculated R^ and K^, and from Eqs. (8) and (9) with 
Koo replaced by max[K d {z)\. In this simulation, R b = 0.2, T b — 5, and u; 0 = 0.8, and the sea 
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surface illumination was taken to be sunny (as defined above). Since R b = R(z — ► z b ), the 
first approach for determining R b is the most straightforward. This gave a smooth profile 
of R b (r) that monotonically approached R b with increasing depth; however, this estimate is 
extremely inaccurate except very close to the bottom, with 41% error at only one optical 
depth above the bottom and 61% error at two optical depths above the bottom. This sharp 
increase in R(z) near the bottom is typical 10 and makes extrapolation of R{z) from mid- 
water depths to the bottom impractical. At all depths off the bottom, far better estimates 
of R b were obtained with Eq. (8) with R^, and K 0 0 calculated from the known water IOP’s. 
The errors at two and one optical depths off the bottom were 11% and 5.8%, respectively. 
Even better estimates of R b at all depths off the bottom, however, were obtained from 
Eq. (8) with R^ determined with Eq. (9) and estimated by max[i7 rf (z)]. This gave an 
error at two and one optical depths off the bottom of 6.2% and 2.8%, respectively. Again, 
extrapolation from mid-depths gave an excellent estimate of R b . 

Values of the bottom albedo from irradiance measurements at a combination of shallow 
and deep-water sites with Eq. (7) were consistently underestimates. Therefore, since 
was generally larger than K d (z), results Were more accurate when the value of K in Eq. (7) 
was determined from a deep value of K d {z ) in the deep-water site, where Kd(z ) ~ AToo, than 
when it was determined from either the deep-water Kd{z) or Kd(z). For example, shown in 
Fig. 2 are R b (r) calculated from Eq. (7) with K replaced by lQ{z ) and by A'oc [determined 
from Kd(r = 15)] for two sunny sky simulations: 1) cu 0 = 0.7, r b = 5, and R b = 0.1, and 2) 
cco = 0.9, r 6 = 5, and R b = 0.2. 

Shown in Fig. 3 are comparisons of R b (r) obtained from E u {z) and E d {z) measurements 
at only the shallow water site [Eqs. (8) and (9)] and i?&(r) obtained from measurements at 
both shallow and deep water sites [Eq. (7)]. The two cases shown are sunny sky simulations 
with r b = 5 and 1) u; 0 = 0.9 and R k = 0.1, and 2) u>q = 0.8 and R b = 0.2. As these examples 
demonstrate, the two-site method was typically far more accurate than the one-site method 
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near the surface, whereas the one-site algorithm usually outperformed the two-site method 
near the bottom. 

Example Rb{r) obtained from L u (z) and E d (z) are shown in Fig. 4 for 1) u ; 0 = 0.9, 
Rb = 0.1, and sunny conditions, and 2) u 0 = 0.9, Rb = 0.2, and overcast conditions. 
Because estimates of Rb with Eq. (11) were typically larger than the actual value, the i?&(r) 
in Fig. 4 obtained with Eq. (11) were calculated with K d (z) rather than the asymptotic K d 
in order to make the estimated values as good as possible. In general, the value of ujq had a 
large impact on the accuracy of both L u -Ed methods, with the best estimates obtained when 
u)q was large. The illumination conditions were also very important for the L u -E d methods, 
with the best results being obtained for overcast conditions. The value of R b , on the other 
hand, had no significant effect on the accuracy of its estimate. Under sunny conditions, the 
one-site method of Eq. (10) and the two-site method of Eq. (11) performed similarly in the 
bottom one optical depth when u 0 was small and in the bottom three optical depths when u> 0 
was large, but Eq. (11) was considerably more accurate than Eq. (10) near the surface. For 
overcast conditions, Eq. (10) performed well at all depths but still underperformed Eq. (11). 

4. Discussion 

Several methods were evaluated here for the determination of R b from common natural light 
measurements. Each method returns accurate values of Rb if implemented very close to 
the bottom. However, because it is difficult in practice to obtain light field measurements 
very close to the bottom, it is necessary to apply these algorithms at one or more optical 
depths off the bottom and, when possible, extrapolate the depth-dependent estimates to the 
bottom. Therefore, it is desirable that the error of the method used is both small and linearly 
decreasing with depth. The methods were all found to be preferable to a straightforward 
extrapolation of R(z) to the bottom, where R(z) = Rb, but differed in their accuracies when 
applied more than one or two optical depths away the bottom. 
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Estimates of Rb can be obtained with Eq. (8) by first estimating /?<*>, either with Eq. (9) 
or by calculating it from known water IOP’s. This method does not require measurements 
at other wavelengths or at another site, and estimates of Rb from this method at one or two 
optical depths off the bottom were generally found to be more accurate than those obtained 
with Eq. (7). Numerical simulations indicated that the use of from Eq. (9) produces 
better estimates of R b than the use of the value of Roo computed from the IOP’s. This is 
because, fortuitously, the error introduced into Eq. (8) by applying it where the light field is 
not well described by the assumed two-mode asymptotic model (see Appendix) is mitigated 
by the deviation in the value of R predicted by Eq. (9) from its true value. Therefore, 
even if the water IOP’s are known, it is preferable to employ Eq. (9), provided the processed 
data has a relatively smooth R b (r) profile. Unfortunately, this method is often inaccurate 
near the sea surface when the bottom signal is not strong. 

If one wishes to estimate Rb from measurements close to the surface and a suitable deep- 
water site is available, then Eq. (7) provides the most reliable method. However, it was fo un d 
that if the deep water site is vertically well mixed, then Eq. (7) should be implemented by 
replacing K^nd. with Koo, which can be directly measured deep in the euphotic zone of the 
deep water site. 

Estimates of R b alternatively can be made from measurements of L u (z) and Ed{z) with 
Eq. (10), provided that the bottom is approximately Lambertian. Equation (10) requires 
that R^ either be analytically computed from local measurements of the water IOP’s or 
measured at a nearby deep-water site. If a suitable second site is readily available, Eq. (11) 
should be used instead since it was found to be generally more accurate and reliable than 
Eq. (10). However, both L u -Ed methods performed well when applied in the bottom half of 
the water column. The inaccuracy of Eq. (10) under sunny conditions makes it unsuitable 
for remote sensing applications, and therefore Eq. (10) offers no advantage over the E u -Ed 
method (which is the more natural approach to in situ R b estimation). Equation (11), on 
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the other hand, shows some promise for remote sensing applications for large u 0 (or for any 
ujq if uj 0 is known). 

Regardless of the method used, the determination of Rb is easiest when the bottom signal 
is strong. Thus, Rb can be obtained most accurately when the water is shallow, the value of 
Rb is large, and the attenuation of the water is low (e.g., over tropical coral reefs or white 
sandy beaches). If the bottom composition is believed to be uniform over a large horizontal 
region, then the determination of Rb should be made at the shallowest depth. 

In practice, the method to use for the estimation of Rb will be dictated by the instrumen- 
tation available and whether or not an appropriate deep-water site exists. Given the choice, 
however, the estimation of Rb should be done with measurements of E u (z) rather than with 
L u (z). The most informative approach would be to use all the methods discussed here and 
intercompare the results. If the estimates agree they can be recorded with great confidence. 
On the other hand, the difference between estimated Rb values from R(z) and R L (z) might 
serve as a crude measure of the degree to which the bottom behaves as a Lambertian surface. 
However, for purposes of modeling the light field, Mobley 22 has demonstrated that it is far 
more important to measure the magnitude of the effective bottom albedo than it is to obtain 
the detailed angular pattern of the bottom bidirectional reflectance function. 
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Appendix A: Derivation of Eq. (10) 

If the optical properties of the water are spatially uniform and there are no internal sources, 
then when cz is at least one optical depth away from a boundary the upward radiance and 
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downward irradiance can be expressed as summations of eigenmodes, 13 

Lu(z) = L(z, -1) = J2 -1) exp i-cz/vj) + C{-Vj -1) exp (cz/vj)], (Al) 

j= 1 
J 

E<t{z) = Y, exp i-cz/uj) + Ci-Vjfai-vj) exp (cz/uj)], (A2) 

3 = 1 

where C(EiZj) are expansion coefficients, are the J eigenvalues of Eq. (1) corresponding 
to the eigenfunctions 23 4>(±Vj, fi), and 

9j(±vi) = [ (A3) 

J 0 

for Legendre polynomial Pj( p). Far from the boundaries, E,i{z) and L u (z) can be ap- 
proximated by retaining only the asymptotic decreasing eigenmode corresponding to the 
largest eigenvalue V\ [i.e., C{uj) — 0 for j > 1 and C(—vj ) = 0 for all j], and since 
R L (z ) = 7r L u (z)/E d (z) the asymptotic value of R L (z ) is R^ = Tx<p{vx, —l)/gx(ui). At 
depths far from the surface but where some influence of the bottom is present, Ed(z) and 
L u (z) are better approximated by also including the eigenmode corresponding to the largest 
negative eigenvalue, —V\, so then 


L u (z ) ~ C(u i)0(zq, — 1) exp(—cz/ui) + C(-Ux)<p(-Ux, -l)exp (cz/ux), (A4) 

E d (z) « C(vx)g\{v\) exp(—cz/vx) + exp(c2yVi). (A5) 


After forming the ratio R L (z) = tx L u (z)/Ed(z), dividing through by 

C{vx)g\{v\) exp(— cz/vx), letting r = C{—u\)/C{vx), and recognizing that 13 R^ = 
<f>(— ux, l)/gx{vx), we find that 


R l {z) 


Rl o + ™r<K~ v - h _1 ) exp {2cz/v l )/g l {vx) 


1 + rR 0 o exp(2 cz/ux) 
Subtraction of R from Eq. (A6) gives, 


(A6) 


nL/ ?1 _ r l = rexp(2cz/^ 1 )[7r0(-t/ 1 , —3 j/giQyi) - 
1 ' 00 l + ri? 0O exp(2c2/n 1 ) 


(A7) 
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A rearrangement of this equation is, 


[R l {z) - R^]exp{-2cz/u x )[l + rR oc exp(2cz/u 1 )] = r[n<f>(-v lt -l)/g l (v 1 ) - R^R^}. (A8) 


Since the RHS of Eq. (A8) is independent of z, the LHS at arbitrary depth z equals that at 
the bottom. Therefore, 


R l (z) = R^ + [R L (z b ) - R^] exp[— 2c(zb - z)/u x ] 


1 + rRp o exp(2 cZb/W) 
1 + rR^ exp(2 czjv\) 


If \rR 00 exp{2cz b /u{)\ C 1 or if (z b — z) is small, then Eq. (A9) reduces to 


(A9) 


R l (z) = + [R f ‘(z t ,) - R^] exp[-2c(z t - z)/v\\. 


(A10) 


which is analogous to our equation for R(z) derived in a similar manner. 17 Rearrangement 
of Eq. (A10) gives Eq. (10). 
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FIGURES 


Fig. 1. Profiles of estimates of the bottom albedo, with a) the estimate equal to the 

irradiance reflectance, R(r), b) Eq. (8) with known i?oo and Koo, and c) Eqs. (8) and (9). The 
simulation was generated for Rb = 0.2, five optical depth water, uiq = 0.8, and 75% collimated light 
at 30° from the zenith. 

Fig. 2. Profiles of estimates of the bottom albedo, i?&(r), with the two-site E u -Ed method of 
Eq. (7) for (left) uq = 0.7 and Rb = 0.2 and for (right) u>o = 0.9 and Rb = 0.1. The simulations 
were for sunny conditions, and the Rb determination was done with K in Eq. (7) replaced by Ka(z) 
(*). and by the deep (asymptotic) value of Kd{z) (-). 

Fig. 3. Profiles of estimates of the bottom albedo, i?b(r), from E u (z) and Ed{z) with the two-site 
method of Eq. (7) (*) and the one site method of Eq. (8) (-). These sunny sky simulations were 
for (left) ojq = 0.9 and Rb = 0.1 and for (right) u o = 0.8 and Rb = 0.2. 

Fig. 4. Profiles of estimates of the bottom albedo, Rb(r), from L u (z) and Ed(z) with the two-site 
method of Eq. (11) (*) and the one site method of Eq. (10) (-). The simulations were for (left) 
sunny conditions, uq = 0.9 and Rb — 0.1 and for (right) overcast conditions, wq = 0.9 and Rb = 0.2. 
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TABLES 

Table 1. Calculations of the bottom albedo with Eqs. (8) and (9) from simulated irradiance mea- 
surements at three, two, and one optical depths above the bottom and extrapolations to the bottom 
of the latter two. 


u>0 

n 

simulation 

Rb 

illumination 

n- 3 

n- 2 

R b (r) 
n- 1 

extrapolation 

0.7 

3 

0.2 

sunny 

N/A 

0.177 

0.189 

0.200 

0.7 

5 

0.1 

sunny 

0.091 

0.091 

0.096 

0.101 

0.7 

5 

0.2 

sunny 

0.170 

0.176 

0.189 

0.202 

0.7 

5 

0.2 

overcast 

0.172 

0.177 

0.190 

0.202 

0.7 
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0.2 

sunny 

0.167 

0.176 

0.189 

0.203 

0.8 
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0.2 

sunny 

0.186 

0.188 

0.194 

0.201 

0.9 

5 

0.1 

sunny 

0.114 

0.106 

0.102 

0.098 

0.9 
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0.2 
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0.207 

0.201 

0.200 

0.199 

0.9 
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0.2 

overcast 

0.207 

0.201 

0.200 

0.199 

0.9 
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0.3 

sunny 

0.296 

0.294 

0.297 

0.300 

0.9 

5 

0.4 

sunny 

0.385 

0.387 

0.395 

0.403 
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Fig. 1. Profiles of estimates of the bottom albedo, i?6(r), with a) the estimate equal to the 
irradiance reflectance, R{r), b) Eq. (8) with known Roo and K 0 c , and c) Eqs. (8) and (9). The 
simulation was generated for Rb — 0.2, five optical depth water, u o = 0.8, and 75% collimated light 


at 30° from the zenith. 




Fig. 2. Profiles of estimates of the bottom albedo, Rb(r), with the two-site E u -Ed method of 
Eq. (7) for (left) ujq = 0.7 and Rb — 0.2 and for (right) ojq = 0.9 and Rb = 0.1. The simulations 
were for sunny conditions, and the Rb determination was done with K in Eq. (7) replaced by Ka(z) 
(*) and by the deep (asymptotic) value of Kd(z) (-). 
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Fig. 3. Profiles of estimates of the bottom albedo, Rb(r), from E u (z) and Ea(z) with the two-site 
method of Eq. (7) (*) and the one site method of Eq. (8) (-). These sunny sky simulations were 
for (left) u>o = 0.9 and Rb = 0.1 and for (right) u>o = 0.8 and Rb = 0.2. 
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Fig. 4. Profiles of estimates of the bottom albedo, from L u (z) and Ed(z) with the two-site 

method of Eq. (11) (*) and the one site method of Eq. (10) (-). The simulations were for (left) 
sunny conditions, cjq = 0.9 and Rb = 0.1 and for (right) overcast conditions, ojq = 0.9 and Rb — 0.2. 
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